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Abstract 

We consider the effects of quantum fluctuations in mean-field quantum spin- 
glass models with pairwise interactions. We examine the nature of the quan- 
tum glass transition at zero temperature in a transverse field. In models 
(such as the random orthogonal model) where the classical phase transition 
is discontinuous an analysis using the static approximation reveals that the 
transition becomes continuous at zero temperature. 
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Spin glasses are models which deserve considerable interest 0. In these systems the 
presence of randomness and frustration can yield very rich behaviour. In particular, there 
is much current interest in the behaviour of glassy systems in the presence of quantum 
fluctuations where the nature of the zero-temperature phase transition is driven by the 
competition between randomness and quantum effects rather than thermal fluctuations 0. 
This makes the order-disorder transition in quantum glasses belong to a new universality 
class. 

Much work has been devoted to the study of mean-field quantum spin glass models. In 
particular, attention has been paid to models with a continuous transition in the Edwards- 
Anderson order parameter. The simplest example in these class of models is the Sherrington- 
Kirkpatrick (SK) model 0] in a transverse field. In this system the critical temperature 
is depressed when the transverse field is switched on and vanishes for a critical value of 
the field Analytical work in the quantum SK model reveals that replica symmetry 

is broken in the quantum glass phase at zero temperature ||. This is an indication that 
quantum fluctuations do not destroy one of the most interesting features in glassy systems, 
that is the coexistence of a large number of phases or states. 

There has been also much recent interest in the study of classical spin-glass models 
with a discontinuous transition in the order parameter. These models are characterised by 
the existence of a dynamical singularity at a temperature above the static transition ||. 
Concerning the statical and dynamical behaviour these models are very good candidates 
for describing real glasses |H|. On the one hand, the statics gives a natural explanation 



for the existence of a thermodynamic ideal glass transition driven by an entropy collapse. 
On the other hand, the dynamics of these mean-field models are described by the mode 
coupling equations introduced to describe relaxational phenomena in glasses |p|-|TTp. In 



mean-field models metastable states have an infinite lifetime, hence dynamics is frozen at 
the dynamical singularity well above the static transition temperature. Below the dynamical 
transition temperature the system gets trapped in states with have larger energy than the 
equilibrium one |12). All these features are absent in models with a continuous transition. 



The purpose of this letter is the study of models with a discontinuous transition in the 
presence of quantum fluctuations at zero temperature. The motivation is twofold. Con- 
cerning the statics we note that the transition cannot be driven in the quantum case by 
an entropy collapse. The reason is that the entropy vanishes everywhere at zero temper- 
ature. Concerning the dynamics we can also expect a quite different behaviour from the 
classical case. In macroscopic quantum systems at T = the dynamics is governed by the 
Schroedinger equation and there is no room for any kind of thermal activated processes. 
It could well be that trapping dynamics in the metastable glassy phase are considerably 
modified or even suppressed in the presence of quantum fluctuation effects. 

The main conclusion of this work is that the glassy scenario presented before is indeed 
suppressed by quantum fluctuations in a certain class of spin-glass models. We will provide 
a general proof for this statement within the static approximation [|J. In models where the 
classical transition is continuous it remains continuous at zero temperature. 

The family of models we are interested in are quantum Ising spin glasses with pairwise 
interactions in the presence of transverse field. These are described by the Hamiltonian, 

i<3 i 

where af, af are the Pauli spin matrices and T is the the transverse field. The indices i,j 
run from 1 to N where iV is the number of sites. The Jy are the couplings taken from an 
ensemble of random symmetric matrices. In the case that the are independent Gaussian 
variables this Hamiltonian reduces to the quantum SK model in a transverse field. If the 
Jij are orthogonal matrices then eq.(|l[) reduces to the random orthogonal model (ROM) 
fT3| in a transverse field. At zero transverse field the models become classical and display 
quite different behaviour. The SK model has a continuous finite temperature transition 
without jump in the Edwards- Anderson order parameter JTJ while the ROM presents a 
strong discontinuous transition where the Edwards- Anderson order parameter jumps to a 



value close to 1 at the transition temperature [13 



In order to solve model ([l|) we apply the Trotter-Suzuki decomposition [14[] and rewrite 



the Hamiltonian in terms of classical spins with an extra imaginary time dimension, 

«.// = A £ J y £ °M + B E + C (2) 

i<j t it 

where the time index t runs from 1 to M and the spins o\ take the values ±1. The 
constants A, B and C are given by A = j-^B = \ \og(coth(^)); C = \og{\sinh{^§-)) . 
Now we apply the replica trick and compute the average over the disorder of the replicated 
partition function, 

Z]= f [dJ] £ expCt Kff) (3) 
HI 0=1 

where f[dJ] means integration over the random ensemble of matrices. This integral can 
be done using known methods in matrix theory |T5|JT3[| . The final result of eq.([jj) can be 
written in terms of a generating function G(x) which depends on the particular ensemble of 
Jij couplings via its spectrum of eigenvalues. For the two examples we will consider in this 



paper we have Gsk(x) = ^ (SK model) and G RO m(x) = \log{^^-A) + |vl + 4x 2 — | 
{ROM model). From (|3|) we get, 

ZJ = J dQdAexp(-NF(Q,A)) (4) 

where 

F(Q, A) = + ^Tr(QA) - i TrG(AQ) - \og(H(A)) (5) 

with Q l a l ' h , A* ^ being the order parameter and the trace Tr is done over the replica and time 
indices. The term H(A) is given by, 

H(A) = £ exp(E E Ki^t + B £ *y a +1 ) (6) 

o ah m W at 

and the free energy is obtained by making the analytic continuation (3f = lim„^ F ^„ A - 
where Q*, A* are solutions of the saddle point equations, A*^ = ^^-(G'(AQ)^ and = 
(a*(T^). The average ((•)) is done over the effective Hamiltonian in (Q). We assume that the 
order parameters Q%, A*^ are independent of the time indices when a but they are only 



translational time invariant if a = b. To study the high-temperature phase and the phase 
boundary of the model we consider a general one step replica symmetry breaking solution. 
We divide the n replicas into n/m boxes K of size m such that m divides n. The saddle 
point solution when a ^ b takes the form = q; A^ b = A if a, b e K and Q l ^ b = A* ^ = 
otherwise. For a = b we take Q^' a = ify_^|,AJ* = A\t-t'\- Going to the frequency space 
we define R p = M" 1 E^o 1 e iuJpt R t , A p = M _1 1 e^*A t and <r p = M" 1 E^o 1 e <w * , *<r t+1 
where cg> p = The free energy (|) is given by, 

s-i M-l i M-l i 

^/ = - a? + E + (™ " -rE G(fiRp) - —G(p(R + (m - 

iV p=0 2 p=0 ^ m 



-G(/3(£ - ?)) + ^(^o) - jT ~ m (x) (7) 



iVi — i 

9(cr p , x) = £ (Ap + MBe-**) \a p \ 2 + {2\)*xa - Aa 2 . (8) 



2m 

where dp(x) = (2n)~2e~~dx is the Gaussian measure and S(x) = E<3- p exp(6(cr p , x)) with 

M-l 

E 

In order to investigate the glassy scenario we compute the static and dynamical transition 
temperatures. Following we expand the free energy ([/]) around m — 1, / = f + (m — 
l)/i + 0((m — l) 2 ) and determine the paramagnetic free energy f and the correction /]_. 
We get, 

Pfo = ~ + E RvK - \ E - J ( A ) (9) 

p p 

(3fi = q\- yG'W + ^ G (^o) - ^G(/3(^o - q) 

+ 1(A) - exp(-/(A)) J dp(x)E{x) log(H(x)) (10) 

where 1(A) = log(H(A = 0)). Note that fo does not depend on q and A as expected for 
the paramagnetic part of the free energy. The static and dynamical transition are obtained 
by solving the saddle point equations ^ = §^ = 0. The static transition appears when the 
free energy / coincides with the paramagnetic free energy / , i.e. f\ = 0. The dynamical 
transition is given by the presence of a soft mode above the static transition and is obtained 
by solving the equation (^)(&) - (0)' = 0. 



The solution to these equations yields the critical temperature and the value of the jump 
of the order parameter q at the transition. In case the dynamical and the static transition 
coincide it can be shown that q = A = and the transition is continuous in the order 
parameter. These three equations are complemented by the saddle point equations for the 
parameters FL,A V i.e. |& = |& = 0. 

We now derive a simple expression for the dynamical transition temperature. The equa- 
tion for the soft mode can be worked out and one finds, 

/Pe-WcrfaRo-qj) jT dp(a;)(<^> - (s ) 2 fz(x) = 1 (11) 

where the average ((•))) is taken over the effective Hamiltonian eq.(|J). For a continuous 
transition (q = A = 0) eq.flTTD can be written in the simple form 

xlG"( X o) = 1 (12) 

where Xo — P{(sq) ~ (^o) 2 ) — PRo is the longitudinal magnetic susceptibility. This equation 
can be solved (for a given G(x)) and yields the critical Xo- I n particular, for the SK model 
G" SK (x) = 1 which yields the result xo — 1 i n agreement with known results In the 
ROM the only solution to that equation is Xo — 00 • F° r a continuous transition this implies 
a divergent susceptibility at the critical field. Using a perturbative expansion in powers of 
1/r it is possible to use eq. fll2"D to obtain T c with reasonable accuracy 0. 

Now we come to the main result of this paper, namely that at zero temperature the 
quantum transition becomes continuous. Then, the static and the dynamical transition 
temperature coincide. We first consider the static approximation where R p = A p = for 
p > 0. Putting R = R , A = A we find the following saddle point equations for R, A, q, A, 

A=f G '( / Jii); R=<< sjnh(^) 

\=t±(G'((3R)-G'(P(R- q )); q =« V °° \\ \')?' ' » (13) 
2 J-oo dp{z)cosh(T) 

where 
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U ° Hodp^aw/iCSoCa:)) 1 J 

with So(x) = (2Ax 2 + /3 2 r 2 )5, T = {b 2 + (3 2 T 2 )^ and b = (2(A - A))*z + (2A)*z. 

Exact expressions are also obtained for the free energies /o, /i and for equation (pp . This 
set of equations can be always numerically solved but explicit results can be analytically 
obtained in the zero temperature limit. Plugging the solution A = u/3, A = v (3 into ([13]) 
and performing the integrals with the saddle point method we find after some lengthy 
computations that u, v and the critical field T c satisfy the equations 

u = -G\ ) (15) 

2 T c -2w v ; 

v = -(G'( )-G'( ] r)) (16) 

It is easy to check that equations (|i"5l) , (|16"[) only admit the trivial solution v = 0. It is also 
possible to show that in case v = also f\ = 0. Because q and R vanish with T and the free 
energy of this solution coincides with the paramagnetic free energy fo we conclude that the 
transition becomes continuous at zero temperature. In order to determine the critical field 
T c we solve eq.([Tl|) in the f3 — > oo limit which yields 

(r c -2 M )- 2 G"( f -^) = l (17) 

with Xo = (r c — 2m) _1 . Equation (|17D together with ( |15|) determine the value of u and 
T c . At the quantum transition point the internal energy is given by U = —T c while the 
entropy is given by S = \G( T * 2n ) — ( r " 2u ) + ^log( r _ 2u )- In case of the SK model 
we obtain u = |,r c = 2 reproducing known results PJT?|. In the case of the ROM we 
obtain u = |,T C = 1. Note that in both models the value of the critical field is given 
by the maximum eigenvalue of the coupling matrix J^. For the ROM, one finds that the 
entropy at zero temperature diverges. This is a consequence of the general failure of the 
static approximation at low temperatures. In figure 1 we show the phase boundaries for the 
dynamical and static transitions in the ROM as a function of the transverse field obtained 
numerically solving equations (0). Both transition temperatures decrease quadratically as 



a function of the transverse field merging into the same point at zero temperature. In figure 



2 we show the Edwards-Anderson order parameter q = < o z > 2 in the ROM as a function 
of r as we move along the static (q$) and dynamical (qo) phase boundaries. 

The static approximation yields inaccurate quantitative results for the thermodynamic 
properties at zero temperature. Nevertheless we expect the order of the transition to be 
correctly predicted. To go beyond the static approximation we should consider all the 
Fourier modes R p , A p in the saddle point equations. This is a non trivial task which remains 
open. 

We stress that the main results of this work are restricted to pairwise interaction models. 



In case of p-spin interaction spin-glass models |T8[ the scenario presented here does not apply 
anymore. For p larger than 2 we have obtained different results indicating that the transition 
does not become continuous at zero temperature. This implies that the quantum transition 
in p-spin interaction models with p larger than 2 belong to a different universality class. 
This is in agreement with known results in case of p-spin models in the limit p — > oo |19| as 



well as in 1/p analytical expansions pO 



Summarising, we have investigated the glassy behaviour in Ising spin glass models with 
pairwise interactions in the presence of a transverse field. In models with a discontinuous 
finite temperature transition we have shown, within the static approximation, that the 
transition becomes continuous at T = and there is no room for a metastable glassy phase. 
This implies that at T = all spin-glass models with pairwise interactions belong to the 
same universality class. 
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Figure Captions 

Fig.l Phase boundaries T S (T) (lower line) and To(r) (upper line) in the ROM in the static 
approximation. At zero transverse field T s ~ 0.0646, To ~ 0.1336. 

Fig. 2 Edwards-Anderson parameter q s (upper line) and qo (lower line) in the ROM on the 
static and dynamical phase boundaries boundaries as a function of the transverse field. 
At zero transverse field q s ~ 0.99983, qr> ~ 0.961. q s and qn vanish linearly with 
at zero temperature. 
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